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Abstract 

Assemblies of circular inclusions with spiraling laminate structure in- 
side them are studied, such as spirals with inner inclusions, spirals with 
shells, assemblies of "wheels" - structures from laminates with radially de- 
pendent volume fractions, complex axisymmetric three-dimensional micro- 
geometries called Connected Hubs and Spiky Balls. The described assem- 
blages model structures met in rock mechanics, biology, etc. The classical 
effective medium theory coupled with hierarchical homogenization is used. 
It is found that fields in spiral assemblages satisfy a coupled system of two 
second order differential equations, rather than a single differential equa- 
tion; a homogeneous external field applied to the assembly is transformed 
into a rotated homogeneous field inside of the inclusions. The effective 
conductivity of the two-dimensional Star assembly is equivalent to that of 
Hashin-Shtrikman coated circles, but the conductivity of analogous three- 
dimensional Spiky Ball is different from the conductivity of coated sphere 
geometry. 

1 Introduction 

Structures with explicitly computable effective properties play a special 
role in the theory of composites. They allow for testing, optimizing, and 
demonstrating of dependences on the structural parameters and material 
properties. These structures are also used for hierarchical modeling of 
more complicated structures and they permit explicitly computing fields 
inside the structure and track their dependence on structural parameters. 
There are several known classes of such structures: the Hashin-Shtrikman 
coated spheres structure [6] and Schulgasser's structures [TT] are probably 



the most investigated geometries of composites. The scheme has been gen- 
eralized to multiscale multi-coated spheres (see the discussion in [HE]), 
coated ellipsoids [3], and the "wheel assembly", studied in JT|. Another 
popular class is laminate structures and derivatives of them, the laminate 
of a rank, which exploit a multiscale scheme in which a course scale lam- 
inate is made from smaller scale laminates in an iterative process. The 
limits of iterations of these structures yields to a differential scheme where 
an infinitesimal layer is added at each step, see 4,8. 

In the present paper, we combine the idea of multi-rank laminates and 
coated spheres, introducing assemblies of circular inclusions with spiraling 
laminate structure inside them. Namely, we study the assemblages of spi- 
rals with inner inclusions, spirals with shells, and assemblies of "wheels" - 
structures from laminates with radially dependent mass fractions. We also 
derive the effective conductivities of complex three-dimensional microge- 
ometries, which we call Connected Hubs and Spiky Balls. The described 
structures model inhomogeneous materials met in rock mechanics, biol- 
ogy, etc. For calculating effective properties, we use the classical effective 
medium theory, see |U [TO] ) coupled with hierarchical homogeniza- 
tion. We study spiral assemblies with inclusions and observe an interest- 
ing phenomenon: a homogeneous external field applied to the assembly 
is transformed into a rotated homogeneous field inside of the inclusions. 
The fields in such structures satisfy a coupled system of two second or- 
der differential equations, rather than a single differential equation that is 
satisfied in Hashin-Shtrikman and Schulgasser's structures. We show that 
the effective conductivity of the two-dimensional Star is equivalent [I] to 
conductivity of Hashin-Shtrikman coated circles, but the conductivity of 
the analogous three-dimensional Spiky Ball is different from the coated 
spheres geometry. 

2 Composite circular inclusion 

Consider an infinite conducting plane with coordinates (x\,X2) and as- 
sume that a unit homogeneous electrical field e = (1, 0) T is applied to the 
plane at infinity, inducing a potential u, 



A structured circular inclusion of unit radius < I is inserted in a 
plane. It consists of a core inner circle of radius ro and an enveloping 
annulus. The inner circle (nucleus) flj = {x : \\x\\ < ro < 1}, is filled with 
an isotropic material of conductivity en. The annulus tt a — {x : ro < 
\\x\\ < 1} is filled with anisotropic material whose conductivity tensor 
S e {r) depends only on radius. The plane outside the inclusion is denoted 
fi„, it is filled with an isotropic material of conductivity a,. 

Effective conductivity of such an assembly is computed by effective 
medium theory. Given the inclusion, we find the conductivity er* so that 
the inclusion is cloaked, 



lim u(x) 



= Xl. 



(1) 



u(x) — xi if | |x| | > 1. 



(2) 
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If this is the case, the outer conductivity <r» is called the effective conduc- 
tivity of the inclusion. The inclusion is not seen by an outside observer, 
therefore the entire plane can be filled with such inclusions, according to 
effective medium theory pi8ll9llT0]. 

We show that the field inside the nucleus Qi has the representation 

u = p (cos(ip)xi + sm(4>)x2) , ||a;||<ro (3) 

where p, ip are constant. An observer inside Qi records a homogeneous 
field similar to the outside field, but rotated by an angle ip. 

The current j, electric potential u, and electric field e in a conducting 
medium are related by equations 

V-i = 0, e = Vu, j = Ke, (4) 

where K is a positively defined symmetric conductivity tensor that rep- 
resents the material's properties. In our assemblage, 

{Oild if x G Qi 
S e if x £Q. a , 
crjd if x G fi* 

Here, Id is the identity matrix. Equations Q are combined as a conven- 
tional second order conductivity equation 

V-A"Vu = 0. (5) 

On the boundaries between different regions, tensor K is discontin- 
uous. But the potential u, the normal current j ■ n, and the tangential 
potential e • t axe continuous, 

[«]t=0, b"-n]t=0, [e-t]t = 0, (6) 

where [.] 1 denotes the jump. For instance, at the boundary of the nucleus, 
we have 

[u]t = lim u(r + e, 0) - u(r - e, 9). 

6->0, C>0 



3 Fields in a spiral assemblage 

In this section we find the fields, currents, and effective conductivity of 
the described assemblage assuming that the eigenvectors of K in Q. a form 
a family of logarithmic spirals. We call the resulting structure the Spiral 
with Core. 

Single Inclusion. Rewrite the problem in polar coordinates r, 9. 
Solving conductivity equation jBJ, we find that the potential in the inner 
and outer isotropic regions is 



\ u(r,9)=Arcos(9) + Brsin(9) if < r < r . . 

\ u(r,9) =rcos(0) if r > 1, {x e 0.) ''' 

A, B are constants. The form of the solution in the outer domain reflects 
the effective medium condition: the inclusion is invisible, and field is 
unperturbed and agrees with condition Q. 
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Figure 1: Spiral with Core 



Assume that the angle <j> of orientation of the principle axes of K in 
fl a to the radius is constant, so that the eigendirections form a family of 
logarithmic spirals. In Q a , K has the form 



K 



K rr K r e 
K r g Keg 



R 



where 

cos(</>) sm((f>) _ in 

— sin(0) cos(<j!>) ' e (72 

<Ti and (72 are positive constant eigenvalues of S e , and <f> is the angle of 
the spiral. The entries of K are: 

K rr — (7i cos 2 ((f)) + (72 sin 2 (f>, 
K r g = (<ti - <r 2 ) cos(^)sin(0) 
ifefl = (7i sin 2 (<?!>) + (7 2 cos 2 (0) 



Equation |5| in the annulus becomes 



( 



id_ 

r dr 



ftrrT; 1" K r g- — 

Or r o9 



+ 



ld_ 
r 89 



d Id 
or r do 



)u(r,0) = O 



(8) 

We separate variables and account for boundary conditions. The solution 
has the form 

u(r, 9) = U{r) cos(f9) + V(r) sin(0). 
and U and V satisfy a system of ordinary differential equations 



LiC/-L 2 V = 0, LiV + L 2 t/ = 
where L\ and L2 are linear differential operators, 



t ( d vr d 
Li = ( — rK r 



1 



v dr 



dr 



Kee 



( d „ d 
V dr dr 



(9) 



(10) 



Notice that equations in Q are respectively the real and complex part of 
a complex- valued differential equation 

(L 1 +iL 2 )(U + iV) = (11) 

Write the current in the spiral as j = J u cos(9) + J v sin(#). We compute 



®L is v - k- ®L 

J u — l\ rr -j- l\r6 , Jv — J\rr T 

dr r dr 



u 



(12) 
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The jump conditions Q at r — ro and r = 1 have the form 

[E7]t = 0, [V]t = 0, [J u ]t = 0, [J„]l = 0. (13) 

Potentials U, V in annulus fi are found by solving 

U = r a (Cicos(/3) - C 2 sin(/3)) + r~ a (C 3 cos(/3) - C 4 sin(/3)), (14) 
V = r a (Ci sin(/3) + C 2 cos(/3)) + r^Gs sin(/3) + C 4 cos(/3)), (15) 

where 

ln(r)Jf rfl 

a = V^iO"2, /O = ^ • 



Determination of Constants. Jump conditions ( 13 \ at the outer 
boundary (r = 1) are 

C 2 + Ca = 0, C 2 - C 4 = (16) 

yielding C 2 = 0, C4 = 0. There remain five unknown constants, k* , Ci, C3, A, 
and B. They appear to be overdetermined by the remaining six bound- 
ary conditions. However, the four boundary conditions ( 13 1 on the inner 
boundary at ro can be reduced to just three conditions, allowing the sys- 
tem to solved. 

The potentials U and V in the spiral material are linked together, 
and knowing either one is sufficient to determine the other. Indeed, upon 
substituting C 2 = 0, C4 = 0, we find 

V(r) = U(r) tan(/3(r)) (17) 

dm = dm tan(/9(r)) + ^Xl+tan^)))^ 
dr dr rK r 



Substituting this into ( 12 1, we see that the same relationship holds for the 
normal components of the current, 

Jv = K r ^V^tan(/3(r)) + ^-K r0 tan 2 (/3(r)) = J u tan(/3(r), (19) 
dr r 



so we have, in the spiral, 

V(r) = Mr)_ 
U(r) J„(r) 



= tan(/3(r)). (20) 



Using (|20f, the current conditions |13j) are redundant, as the two con- 
ditions on the current are simultaneously satisfied. The four boundary 
conditions can be rewritten as three, 

[U}±=0, [V}±=0, tan(/3(r )) = ~. (21) 

With this observation, the linear system can be solved for the remain- 
ing unknowns. Among them is the effective conductivity of the outer 
region, a, , which is treated as an unknown. Notice that the annulus with 
spiraling material does not perturb the outside field, and contains a ho- 
mogeneous inner field that is directed in a different direction than the 
outer homogeneous field. 
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The value of a* is obtained by solving the system is the effective con- 
ductivity of the structure; it is given by an explicit formula 

gggifc? - 1) + a t ^ai(J2{rl + 1) 

- 1) + y/wHrZ + 1) ' 

where 

2a/(Ti(J2 

7 = • 

sin(0) 2 <Ti — sin(0) 2 (72 — <ti 



Figure 2: Spiral Assemblage 

Effective Medium Theory. The Spiral with Core structure, as 
described above, extends from the origin to a radius of one. But this 
radius is arbitrary in an infinite plane. The spiral can be scaled up or 
down with respect to the radius and will still solve the same problem. As 
the spiral leaves the outside field unperturbed, placing multiple spirals side 
by side will result in every object placed being rendered invisible to the 
field, with each having an inner homogeneous field directed in a different 
direction than the applied field. While none of the myriad inclusions are 
detectable by perturbations of the outside field, the field is differently 
directed almost everywhere. 



4 Extreme Spiral Geometries 

Extremal Spiral Angle. By design, the core of the spiral object contains 
a uniform electric field directed in a different direction than the outside 
field. The potential in inner region (nucleus) has the representation 

u(r,8) = Arcos(6») + Brsm(9). 

The field inside will make an angle of T = tan -1 (B /A) with the uniform 
field outside the spiral. A and B depend on the spiral angle (f), so a 



G 



substitution shows T = 7|<72 — where 

ln(ro) cos(c6) sin(c6) 



7 = 



J (c4)cn - sin (0)cr 2 - en 



(23) 



Let us find the angle cj> which maximizes the rotation T. A straight- 




Figure 3: Changing Direction of the Electric Field 

forward calculation shows that for a given en and 02, T is maximized by 
choosing 

4>o = arctan(y^CTi7^2~)- (24) 
The maximal angle T max for a given en, CT2, and ro is 



T„ 



■ ln(r ) 



11 

a 2 



(25) 



Here, ro < 1, the radius of inclusion. It's clear that the more anisotropic 
the spiral and the smaller the inner radius ro are, the larger the resulting 
twist inside the spiral's core is. 



Laminates. The conductivities en, en in (251 describe the conductiv- 



ities of the outer spiral material in the Spiral with Core structure. If this 
anisotropic material is a laminate made of materials with conductivities 
ki,kz and volume fractions mi, 7712, then en, CT2 can be written as the 
geometric and arithmetic means of the conductivities, [4] 



en = mifei + 7712^2, cr 2 



(26) 



/ mi mi • 
V fci fc 2 

For given materials with conductivities ki, fe, we can optimize T with 
respect to mi, 7712. The maximal angle of rotation is obtained for 



mi 



m,2 = 



(27) 



It is equal to 



T 

J- r. 



-lln(ro) 

4 (fci +fc 2 )\/fcife 



(28) 
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In particular, T = n (and the current inside goes in the opposite direction) 
if 

ln M = - 47r (fci _ fca)2 . (29) 

For instance, if ki = 1 and = 100, the inner field will be directed 
opposite to the external field when ro = 0.274. 

Values of T greater than 27r are also possible. The relative angle 
between the incident current and the inner current is the value of T taken 
mod 2n. 



5 Derivative Assemblages 

The parameters in the Spiral with Core structure can be modified to 
obtain effective conductivity of similar conducting assemblages. 

Hashin-Shtrikman Coated Circles. The classical example of the 
Hashin-Shtrikman geometry[3] is an isotropic circle surrounded by an 
isotropic annulus 



Figure 4: Left: Coated Circles. Right: Schulgasser's Geometry (orange slice) 

It is obtained from the Spiral with Core by setting o\ — <T2. The 
outer spiral layer becomes an isotropic shell, and the effective conductivity 
coincides with Hashin-Shtrikman's result 

, (gj +0-1) +rf ) {a l -cri) 

kh„ = <n-, : c Yi V- ( 30 ) 

Notice that m = tq is the volume fraction of the core material. 

Schulgasser Structure. Schulgasser suggested another clas- 
sical symmetric geometry . It is a radial laminate of two materials. 
Schulgasser's structure is obtained from the Spiral with Core by setting 
ro = 0, (j> = 0. The inner isotropic circle disappears and the logarith- 
mic spiral laminate is straightened into a radial laminate. The effective 
conductivity agrees with |11| . 

ksch = V "!^- (31) 

Orange with Core. This object has an inner circle of isotropic 
material surrounded by an annulus made up of a radial laminate. 

The Orange with Core is obtained from the Spiral with Core by set- 
ting cj> — 0. The logarithmic spiral laminate is straightened into a radial 
laminate. The effective conductivity of the Orange with Core is 
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Figure 5: Left Orange with Core. Right: Orange with Shell 



ffiK (l + r'i K ) + <t 2 (l ■ 



K c = (Jl ; =— ; ?r^> K=\ • (32) 

o-ik (1 + rg K ) + <7i (1 - rg K ) V' 71 
Note that this formula can be written in terms of volume fractions by 

observing that the volume fraction rrii of the inner material is rrii = Vq. 
Orange with Shell. This object contains an inner radial laminate 

surrounded by an isotropic material. Denote the conductivity of the outer 

isotropic shell by ai, and the conductivity tensor (in polar coordinates) 

for the inner laminate by 



K = 



ay 
ag 



The effective conductivity can be found by a straight calculation. It is 



k a = 01 



y/Wr{r 2 + 1) + o-i(l - rl) 



(33) 



' \A 7 »°>(1 - r o) + ff i ( r o + !) ' 
Indeed, the material in an inner circle is a Schulgasser radial laminate, 
which has effective property k„ = yJo r og. Thus, this inner material can 
be treated as an isotropic material with conductivity y/o r ae. The effective 
property of the Orange with Shell can be obtained by substituting ai = 
^a r ag into the equation for the effective property of Hashin-Shtrikman 
coated circles. 

Basic Spiral. This object is simply the spiral material centered 
around the origin. 





Figure 6: Left: Basic Spiral. Right: Spiral with Shell 

The Basic Spiral can be obtained from the Spiral with Core by setting 
ro = 0. The effective conductivity of the spiral is 

k sp = \Jo\(J2- (34) 
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As one might expect, it coincides with the effective conductivity of Schul- 
gasser's geometry. 

Spirals with Shell and Core. The Spiral with Shell contains a circle 
of spiral laminate surrounded by an isotropic shell. As we mentioned, no 
more work is needed to be done to calculate the effective property of it. 
Since the effective property of the Basic Spiral is the same as the effective 
property of Schulgasser's Structure, the effective property of the Spiral 
with Shell is the same as the effective property of the Orange with Shell, 
k s - ( |33[ ) For completeness, we also mention here again the Spiral with 
Core discussed in Section 3 that is an isotropic circle surrounded by a 
spiral laminate. It is a field rotator, and its effective property is given by 
§2). 

Insulated Geometries. Consider the asymptotic case: assume that 
the spiral is a laminate made from two isotropic materials, and the con- 
ductivity of one material approaches zero. 

If one of the conducting materials is replaced by an insulator, the 
conductivities of both Schulgasser structure and the spiral are zero. The 
spiral with shell structure can be insulated by replacing one of the ma- 
terials in the spiral laminate with an insulator. The resulting structure 
should behave as an annulus; the inner spiral acts as an isotropic insula- 
tor, and all current should passes through the isotropic shell. Indeed, the 
resulting effective property is 

, o v ; 1 — r o 

k hs = hm k hs = — - — jo-i. 
<r 9 ->o 1 + 

This structure has the effective property of an annulus with insulated 
nucleus, as expected. 

If the conductivity of one of the materials in the spiral laminate is 



zero, the spiral with core structure has effective property (see ( 22 l ) 



fc2=Iimft.- ^cosfo) 



ov cos 2 (</>) — Oi ln(ro) ' 



The effective property of if the insulated orange with a core obtained from 
this one by setting = 0, 



ai ln(r ) 



Wheel. Insulators can be used in laminates to control the variation of 
the current material in the layers. Consider the following Wheel geometry. 
A central core is surrounded by an annulus that consists of conducting 
radial spikes surrounded by an insulator. The thickness of an individual 
spike is constant, and the volume fraction of conducting material in the 
annulus decreases with radius, because the circumference length linearly 
increases. The conductivity tensor is 



K = 



o-i/r 




(35) 



A spike conducts in the radial direction but does not conduct in the cir- 
cumferential direction. Call this material a spoke material. To relate cri 
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to the conductivity a of the spoke material, we can write afiro = o\. The 
parameter /i £ [0, 1] shows the relative thickness of the spikes at the radius 
ro- 

This time we are dealing with a material with anisotropic conductivity 
that varies with radius. The solution to the potential inside the spoke 
material satisfies 151. The current is constant inside each spoke, therefore 
the potential is 

u = {Ar + B) cos(0) + (Or + D) sin(0). (36) 

The wheel inclusion has an inner isotropic core and an outer isotropic 
shell, connected by spokes. It conductivity depends on radius as follows 



Figure 7: Left: Wheel. Right: Star 



K = 




if 

if 

if 
if 



r < r 

ro < r < n 

ri < r < 1 
r > 1 



(37) 



Since the potentials in each region are known, a procedure similar to 
that in Section 3 gives that the effective conductivity of the Wheel: 



where 



A 



k m heel = 0"2 



A + B 
B-A' 



(38) 



B 



r x (71(72 

Ti<Ji(J2 + 02fl + (7l<7; — r tJ 2 Vi 



When ri — > 1, the outer conducting layer disappears and we come to 
the structure that we call Star. Its effective conductivity is 



0101 



(39) 



01 + <7i(l — r ) ' 

Star assemblages can model structures of hubs connected by conducting 
strands in an insulating space. If the conductor in inner circles and spoke 
is the same, &i = a, we write the effective property a s tar in terms of the 
volume fraction m of conducting material a, 



k s tar{ro,m) = 



am(r a ) 
2 - m(r ) 1 



where m(ro) = r + 2/iro(l — ro). (40) 
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If p = |, then ro = m and fcstar coincides with the effective con- 
ductivity of coated circles, therefore it is an optimal structure for the 
resistivity minimization problem [I] along with the coated circles [B]. The 
value n=\ makes intuitive sense. One can see that if the spokes in the 
Star cover only half the circumference of the inner sphere at ro, then the 
current density in the isotropic region of the star will be half that in the 
spokes region of the star. If two orthogonal currents are separately applied 
to the Star structure, the sum of their energy density is constant through- 
out the structure. This satisfies a necessary and sufficient condition for 
minimization of resistivity [4] 



6 3D radially symmetric assemblages 

The same approach can be applied to three-dimensional radially symmet- 
ric inclusions. Let us consider spherical inclusions, embedded in an infinite 
three-dimensional isotropic conducting material with undetermined con- 
ductivity. A uniform electric field is applied to the plane. The equation 
(JsJ is used to solve for the potential in each material. The boundary con- 
ditions |6]| are similar to the two-dimensional case. Applying a uniform 
current along the z-axis results in the electric potential 

u = pcos(0) (41) 

outside the inclusion. This solution is obtained by separation of variables 
in equation (|5| in spherical coordinates. 

Hub. We construct an analogue of the two-dimensional spoke mate- 
rial, which consists of radial spokes of constant thickness. The relative 
proportion of conducting material is inversely proportional to the square 
of the radius, and the tensor in spherical coordinates is 



p 2 



1 





(42) 



We choose p £ [0, 1] so that the term ppo is the surface area covered by 
the spikes at the inner radius po- 

The solution for the potential inside this material is 

u= (C 1 + C 2 p)cos(cj>). (43) 

Using this material, we construct a Hub structure. The Hub is the three- 
dimensional analog of the two-dimensional Star structure. It model in- 
clusions spread through an insulating space and connected by conducting 
strands. The Hub consists of an inner isotropic core of conductivity Oi 
out to a radius of p, enveloped by a spoke material. 

The effective conductivity of the Hub is computed similarly to that of 
Star. It depends on two structural parameters, po and p and is equal to 

a t ap%p . , 

Khub = —7- r— ■ (44) 

<Ji(l — po) + crppo 

The volume fraction of the Hub is 
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Figure 8: Left: Hub. Right: Spiky Ball. 



m = rg + 3/irg(l-ro). (45) 

If, by analogy with the star, we take /u = | , this relation gives r% — m, 
and the effective conductivity is 

k hub = a ™ ' . (46) 
3 — 2ym 

Unlike the two-dimensional case, this conductivity is different from the 
conductivity of coated spheres. 

Spiky Ball. The Spiky Ball structure is a generalized hub, in which 
the spikes have variable cross-section. The tensor of this generalized ma- 
terial is 



K = °Vfo 
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n" 

0. 

The effective conductivity of this assemblage is 

aiapofj,(n — 1) 
Po ) + <Wo ( n - !) 



(47) 



k *b = — -^=TT -■ ( 48 ) 



Notice that the Spiky Ball, Hub, and Star require that the current 
passes sequentially through insulated spokes and a homogeneous central 
circle (sphere). Because of this, their effective conductivity are of har- 
monic mean type, 

lab , 
— = - + —, (49) 
fc* a Oi 

where a and b are constants depending only on the geometry. 



7 Conclusion 

In constructing the Spiral with Core structure, we constructed a field 
rotator, e.g. a conducting structure that rotates an incident uniform elec- 
tric field while remaining undetectable. Optimization results show how to 
maximize the angle of rotation induced by this field rotator. By modifying 
the Spiral with Core, a large number of other conducting structures can be 
constructed and their effective conductivities computed. The analogous 
3D Hub and 2D Star structures behave differently in the sense that the 



13 



2D star corresponds to minimal effective conductivity while the 3D Hub 
does not. 
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